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ENHANCED ELLIPTIC GRID GENERATION 

CROSS REFERENCE TO RELATED 
APPLICATIONS 

5 

This application claims the benefit of U.S. Provisional 
Application 60/425,750 filed Nov. 7, 2002. 

ORIGIN OF THE INVENTION 

to 

The invention described herein was made by an employee 
of the United States Government and may be manufactured 
and used by or for the Government for governmental pur- 
poses without the payment of any royalties thereon or 
therefor. 15 

TECHNICAL FIELD 

The present invention is a method for eliminating require- 
ments for parameter inputs for generalized grid generation in 20 
modeling of engineering systems. 

BACKGROUND OF THE INVENTION 

A large amount of effort has been devoted to developing, 25 
enhancing and using grid generation techniques, through 
solution of elliptic partial differential equations (“PDEs”). 
Elliptic grid generation methods generally focus on devel- 
oping body-conforming grids around bodies for simulations 
of external fluid flow. The grids thus generated are smooth, 30 
having at least continuous first and second derivatives, 
appropriately stretched or clustered, and are orthogonal over 
most of the grid domain. Inclusion of inhomogeneous terms 
in the PDEs allows a grid to satisfy clustering and orthogo- 
nality properties in the vicinity of specific surfaces in three 35 
dimensions and in the vicinity of specific lines in two 
dimensions. 

Following the work of Thompson, Thames and Mastin, 
Jour. Computational Physics, vol. 24, 1977, pp. 274-302, 
three-dimensional governing equations for elliptic grid gen- 40 
eration are often expressed as: 

(1A) 

n«+n yy +n 22 =2(5. i i^)=-«i««(’i- i i,)“p{-'i,tn-'n,o. ( JB ) 

45 

^+C w ,+^=*^4,0=-AS«»(C-?,)ex pl-^,1}, (1C) 

where i;, tj, and £, are generalized curvilinear coordinates, x, 
y and z are Cartesian coordinates, and P(1;,T|,Q, Q(£,Tj,Q and 
R(£,,r|,Q, are inhomogeneous terms, a,, b i( c„ d„ e, and f f are 50 
manually selected constants, and the subscript “i” refers to 
a particular boundary component associated with the prob- 
lem. 

In a two dimensional study by Steger and Sorensen, Jour. 
Computational Physics, vol. 33, 1979, pp 405-410, the 55 
authors use the following governing equations, 

^+4 w =-n i -jgn(il-ri,)exp{-</,tn-i1.1}. (1°) 

n«+ 1 i B .=- c .-xs"(’i-nf)exp{-</,tii-n,fl. (iE) 

60 

for a given T) boundary. The quantities a,- and c,- are gener- 
alized to functions a,-© and c f (i;), respectively, and the 
values of these quantities are computed as part of the 
solution by requiring a specified spacing between a given q 
boundary and an adjacent grid line, and grid orthogonality at 65 
this q boundary. In any two dimensional problem, the decay 
parameters, such as d f , must be prescribed or manually 


2 

inserted for each of the boundaries in any coordinate direc- 
tion. However, as will be seen in the subsequent develop- 
ment, these values of d,- are coupled with the values com- 
puted for the quantities a, and c, respectively so that explicit 
prescriptions of values for the parameters d,- are in conflict 
with the values computed for a,- and for c,. Further, the 
process of selecting the two values for the parameters d, for 
two opposing boundaries, and, by extension, four values for 
four boundaries in a two dimensional problem, is cumber- 
some for static grids and is infeasible where dynamically 
changing grids are required. In three dimensions, the param- 
eter values need to be prescribed for six boundaries, one for 
each of six boundaries. 

What is needed is an approach that provides an automatic 
procedure for generating an elliptic grid that does not require 
manual insertion or user prescription of these decay param- 
eters for a two dimensional or three dimensional grid 
generation problem. Preferably, these decay parameters 
should allow for a variable rate of decay from different 
points on any grid boundary, should arise automatically in 
the formulation and solution of the problem and should 
permit an interpretation in terms of one or more physical 
quantities associated with the problem. 

SUMMARY OF THE INVENTION 

These needs are met by the invention, which provides a 
process for generating an elliptic grid in which ( 1 ) grid 
points tend to cluster near a boundary at a desired rate, which 
may vary from point to point where one or more coordinate 
variables may undergo a relatively large change in value), 
( 2 ) grid lines corresponding to a constant value of a coor- 
dinate are approximately parallel to or perpendicular to the 
local boundary line, and (3) user prescription or manual 
insertion of parameters to achieve a deshed grid behavior in 
terms of clustering and orthogonality near boundaries is not 
required (or allowed). The process includes the following 
steps: 

providing defining equations, valid near at least one 
boundary segment in a generalized coordinate system, of a 
selected grid system, where each of the defining equations 
has at least two independent Cartesian coordinate variables, 
has at least one generalized coordinate as a dependent 
variable, and comprises a partial differential equation, 
expressed in at least one generalized coordinate; 

providing a selected group of boundary constraints for the 
grid system, valid near the at least one boundary segment, 
where a decay parameter for at least one of the generalized 
coordinate dependent variables near the at least one bound- 
ary segment is determined as part of a solution for the grid 
system, rather than being prescribed initially; 

providing defining equations and selected boundary con- 
ditions, having at least two independent coordinate variables 
and at least one dependent variable, for steady state heat 
transfer on a long thin fin, and providing a correspondence 
between the at least two independent coordinate variables 
for the grid system near the at least one grid boundary 
segment with the at least two independent coordinate vari- 
ables for the heat transfer problem; 

providing a correspondence between a selected power of 
at least one heat transfer coefficient for the heat transfer 
problem and at least one decay parameter for the grid system 
near the at least one grid boundary segment; and 

determining a solution of the grid system near the at least 
one grid boundary segment that incorporates at least one 
boundary constraint comprising the at least one decay 
parameter determined for the grid system. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a flow chart of a procedure for practicing the 
invention. 

FIGS. 2 and 3 illustrate grids, computed using the inven- 
tion, for a two dimensional annular region and for a two 
dimensional convex region, illustrating clustering of grid 
points near an inner boundary and near an upper boundary, 
respectively. 

DESCRIPTION OF BEST MODES OF THE 
INVENTION 

The equations (1A), (IB) and (1C) are modified here in 
the context of equations (ID) and ( 1 E) and are written in the 
following form, as a six-equation set for each of the (^,r),Q 
boundaries, where a t> =a fci (ti,C)c fcf =c fc ,<^,Q, and e* =e M (i 
r|) (k=l, 2, 3). The decay parameters b,-, d, and f, are 
(positive) constants for any given boundary segment and are 
expressed as parameter functions, b,(t|,Q, d|(£j,Q, and f,(^, 
ri). Without loss of generality, one can assume that near a 
given ^-boundary segment i, ^-^,-,§0, in a selected region 
on one side of this boundary segment, where sgn(^-^,)>0. 
Treatment of a situation with sgn(i;-i;,)<0 is analogous. In a 
region close to this boundary segment, where b,(T|,Qli;- 
iy«l, the defining equations and the nonhomogeneous 


terms have the forms 

(2A) 

n„+n w h = =?ift,n,0. (2B) 

(2C) 

(2D) 

(2E) 

9i(4.n,D=-c,.i(n.?)sg"(^,« x p{-6i( T i,0^-^fil, (2F) 

=-c 1 ,(n.?Rc 1 ,(n,D(4-^,), (2G) 

(2H) 

«-e 1 , i (T 1 ,Q+ eu <Tl,0(^,.), (21) 

In a region close to an q -boundary segment, where d f (£,,Q 
lq-q,l«l and q-T|,S0, the defining equations and the 
nonhomogeneous terms have the forms 

^,yy^ Z =P2^A,Q, OA) 

n«+Tl n .+h 5i =92(4.h.O, (3B) 

^„+^r 2 (^n,0, (3C) 

P2(^’1.0=-O2.i(^Dsgn(h-hf)expW,<t0lri-il.l}. (3D) 

(3E) 

92(^'n,C)=-C2.,(^CB«"( T l-'ni)exp{-«/,(^,Qlri-Ti,J}, (3D 

=-c 2-i (^0rf,<5,0(n-Ti,), (3G) 

r 2(^'n,0=-«2.i(5.0x«"( T i-ni)e x p{-<i,(5,0tn-Ti,i}, (3H) 

^e 2 /^ 2 /^)d,(i0(nn,l (3i) 


In a region close to an ^-boundary segment, where f,(^,q) 
l£-£,-«l and £-£,-§0, the defining equations and the non- 
homogeneous terms have the forms 




(4A) 


n«+n»'n«=?3(in.O. 

(4B) 

5 

C„+C w +^='-3(tn.O. 

(4C) 


f3(^ri.Q=- a 3.i(tilh'«n(GC,-)exp{-/ / ft,ri)^-< i l}, 

(4D) 


“-%,i(in)+% t ,Cn)/^n)(GQ, 

(4E) 

10 

?3(5.T|.0=-%.f(^’l)s*n(Wi)exp{-M,Tl)C-? 1 l}, 

(4F) 


=-c 3i ,(^nFc3, i (^ii)/i(5,ii)(C-?i), 

(4G) 


r3(5,n,0=-e 3 .,<til)^"(C-C,) ex P{-/.(tri)l?-?.l}, 

(4H) 

15 

=-e 3i ,(tn)+e3.,<iill/,(iil)(C-?,), 

(41) 


Where the preceding approximations, near a boundary seg- 
ment C,=C,j, for example, are used, the defining equations for 
20 ^ become 

L»+4 w +4a-“3.,(^n)M. 1 iX^i)=-fl3,,(^it>g«(C- 

Q, (5A) 

n«+ny y +hH-c3,,(4DM.'n)K-?,-)=-c3,,€,n)sgn(C-- 
25 Q, (5B) 

(5C) 

If one ignores the nonhomogeneous terms on the right hand 
side of Eq. (5C) and ignores the dependence upon z (or Q, 
30 this relation is seen to represent a steady state heat transfer 
equation for a long, thin fin of width or height 2L, with 
corresponding “temperature” 0=UCi, discussed, for 
example, by V. S. Arpaci in Conduction Heat Transfer, 
Addition Wesley, Reading, Mass., 1966, pp. 145-147 and 
35 201-205: 



(6) 

m 2 = 2 h/k8. 

(7) 


40 

where h is a heat transfer coefficient in a selected (z-) 
direction, k is a thermal conductivity coefficient and 8(«L) 
is thickness of the fin. Similar equations apply for 0=U^i, 
and 0 = 11 - 11 !. 

45 The heat transfer coefficient h corresponds to or is pro- 
portional to, for example, a decay parameter, such as the 
coefficient e 3 ,(i;,T|)-f;(i;,Q in Eq. (41). Similar equations are 
developed for the choices 0=U^ii or 0=T|-T| j . 

Where the heat transfer coefficient h in the z-direction is 
50 small, the thermal gradient is correspondingly large normal 
to the corresponding boundary in the xy-space, which 
requires a close clustering of isothermal lines, or of the 
corresponding time evolving grid lines. This provides a 
physical basis for the observation that clustering near a 
55 given boundary increases as a decay parameter (e.g., m 2 in 
Eq. (7)) decreases, and inversely. However, for a grid whose 
grid lines show a general dependence upon both coordinates, 
x and y (as in FIGS. 2 and 3, discussed in the following), a 
similar conclusion is arrived at using the following physical 
60 argument. As the heat transfer coefficient h in the z-direction 
decreases, the temperature gradient near a y-boundary 
increases correspondingly, for a given heat flux, and this 
leads to a denser clustering of isothermal lines (or, equiva- 
lently, of grid lines for the time evolving grid). 

65 Boundary constraints, valid in a one-sided neighborhood 
of each of the four (two dimensional) or six (three dimen- 
sional) boundary segments, are incorporated by applying 
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Green’s theorem in three dimensions for each of the Eqs. 
(2E), (3G) and (41) for the q and boundaries, respec- 
tively: 

J, 5 0 e/dn) 4 a=J v {-fl 1 ,(11,06,(0,0 

0}4r, (8) 

l s @Bldn)da=S^-c 2 /$,Qsgn(T\~T\ i )+c 2 /k'Qd i (l,&) 

0}4t, (9) 

; s (a0/an)30/3H)4<j=] v {-e3 1 (^,ri)ign(C-</)+e3,,(^,Tl)/, 

ft,n)0}}4t, (10) 

where n refers to a direction that is locally normal to a 

surface S representing a totality of six surfaces including the 
boundary segments of interest and V is a volume enclosed 
or defined by the totality of these six surfaces. These 
integral-type boundary constraints can be used to calculate 
the decay parameter analogs, a M (r),Qb,-(r|,Q, c 2 ,(r|,Qd,(i;,Q 
and e 3 ,(T|,i;)f i (i;,T|). When expressed in terms of the gener- 
alized coordinates (£,T|,Q, the boundary constraints set forth 
in Eqs. (8), (9) and (10) are transformed, for C, for example, 
as follows. 

j s I do=jJdQ/t}n)da=js(d{Jdii)dc ( 11 ) 

The integral j s I do in Eq. (14) can be written as an 
algebraic sum of six integrals, evaluated over the indicated 
boundary segments: 

U do-tJ^ ma J da-j^ ni J 

do-J^Jdo-J^-Jda, (12) 

where the surface configurations cm ax, qmin, etc. represent 
the corresponding boundary segments that together make up 
the surface S. For the first and fourth integral pair, the second 
and fifth integral pair, and the third and sixth integral pair in 
Eq. (15), the following respective relations may be verified: 


J 5 da=Jf(l/7 v /a ll )a 1 3[(Ti2+y n 2 +^ 2 )(x5 2 +>’ t 2 + 
z™ 2 )] 1/2 4(4n. 

(13) 

V da= L(l/7/a 22 )a 23 [(r 5 2 +z n 2 )(.v c 2 +y t 2 + 
z^d^ 

(14) 

J ? / da=Js(llJ)[a 33 (x^ 2 +y^ 2 +z^ 2 )(x i 2 +yi L 2 +zZ‘)] 11 
2dr\d^, 

(15) 


(16A) 

fI 22=T(ll, 2 tr)v 2|T lz 2 ), 

(16B) 

a33=AC, 2 +C/+C/), 

(16C) 


(16D) 


(16E) 

«23= 22 bht+nA +T i 2 Q. 

(16F) 


where J=J((x,y,z)/(^,q,Q) is a Jacobian of the transformation 
(x,y,z)- ->(^,i|/,Q. 

Equations (13)— (15) can be used to express the boundary 
constraints in the computational space (generalized vari- 
ables (!;,\|/,Q), analogous to the Eqs. (8)— (10) expressed in 
Cartesian coordinate space. The defining equations in com- 
putational space that are solved, subject to the boundary 
constraints set forth in Eq. (41) for a ^-boundary, are 

«n^,«l^2X,-. nn ^3-r..);cl 2 {“i2Xi,5T,+“i3X,,5c MX 23Xi, 

(17) 

x~x, y or z, (18) 


6 

FIG. 1 is a flow chart of a suitable procedure for practicing 
the invention for a two dimensional or three dimensional 
grid system. In step 11, the system provides defining equa- 
tions, valid near one or more grid boundary segments in a 
5 generalized coordinate system, of a selected grid system, 
where each of the defining equations has at least two 
independent Cartesian coordinate variables, has at one gen- 
eralized coordinate as a dependent variable, and comprises 
a partial differential equation, expressed in at least one 
to generalized coordinate. 

In step 13, the system provides a selected group of 
boundary constraints for the grid system, valid near the one 
or more boundary segments, where a decay parameter for at 
least one of the generalized coordinate dependent variables 
15 near the one or more boundary segments is determined as 
part of a solution of the defining equations, rather than being 
prescribed initially. 

In step 15, the system provides defining equations and 
selected boundary conditions, having at least two indepen- 
20 dent coordinate variables, for steady state heat transfer on a 
long thin fin, and the system provides a correspondence 
between the at least two independent coordinate variables 
for the grid system near at least one grid boundary segment 
with the at least two independent coordinate variables for the 
25 heat transfer problem. 

In step 17, the system provides a correspondence a 
selected power of at least one heat transfer coefficient for the 
heat transfer problem and at least one decay parameter for 
the grid system near the at least one grid boundary segment. 

In step 19, the system determines a solution of the grid 
system near the at least one grid boundary segment that 
incorporates at least one boundary constraint comprising the 
at least one decay parameter determined for the grid system. 
FIG. 2 illustrates a result of application of the invention 
35 to a two dimensional annular region to provide a grid in 
which grid points cluster near an inner boundary of the 
annulus. Where a geometrical system, such as an annulus, 
evolves with time, a grid according to the invention is 
developed at each of a selected sequence of times, with the 
40 parameters subject to the boundary constraints being 
allowed to vary from one time to another time. Each of these 
grids can be used to perform a finite element or finite 
difference analysis on the geometrical object that represents 
the time evolving system at one of these times. 

45 FIG. 3 illustrates a result of application of the invention 
to a two dimensional convex-concave geometry region to 
provide a grid in which grid points cluster near an upper 
boundary of the region. 

A decay parameter, such as e 3 ,(^,ri)-f,(^,q), may vary 
50 with one or more of the generalized coordinates, such as E, 
and/or r\, rather than being constant; and this variation is 
determined as part of the solution of the grid problem, rather 
than being prescribed initially by the user. A grid solution 
can be determined for a temporally constant environment. 
55 Alternatively, a time evolving environment can be allowed 
to vary at each of a sequence of times, and a grid solution 
can be determined for each of this sequence of times, using 
the preceding analysis at each of these times. 

The preceding analysis has focused on neighborhoods of 
60 the grid boundary segments. As noted in the preceding, in an 
interior region, far from the grid boundary segments, the 
defining partial differential equations (PDEs) become homo- 
geneous, and standard analysis of elliptic PDEs is applied to 
determine an interior solution, which is automatically 
65 matched in the solution process across a selected interior 
boundary to the solution obtained for the grid boundary 
segments. 
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What is claimed is: 

1. A method for generating an elliptic grid of coordinates 
in two dimensions or in three dimensions, the method 
comprising: 

providing defining equations, valid near at least one 
boundary segment in a generalized coordinate system, 
of a selected grid system, where each of the defining 
equations has at least two independent Cartesian coor- 
dinate variables, has at least one generalized coordinate 
as a dependent variable, and comprises a partial dif- 
ferential equation, expressed in at least one generalized 
coordinate; 

providing a selected group of boundary constraints for the 
grid system, valid near the at least one boundary 
segment, where a decay parameter for at least one of 
the generalized coordinate dependent variables near the 
at least one boundary segment is determined as part of 
a solution for the grid system, rather than being pre- 
scribed initially; 

providing defining equations and selected boundary con- 
ditions, having at least two independent coordinate 
variables and at least one dependent variable, for steady 
state heat transfer on a long thin fin, and providing a 
correspondence between the at least two independent 
coordinate variables for the grid system near the at least 
one grid boundary segment with the at least two 
independent coordinate variables for the heat transfer 
problem; 

providing a correspondence between a selected power of 
at least one heat transfer coefficient for the heat transfer 
problem and at least one decay parameter for the grid 
system near the at least one grid boundary segment; 

determining a solution of the grid system near the at least 
one grid boundary segment that incorporates at least 
one boundary constraint comprising the at least one 
decay parameter determined for the grid system; and 
displaying the solution of the grid system. 

2. The method of claim 1, further comprising selecting 
said grid system decay parameter to be directly proportional 
to said heat transfer coefficient. 

3. The method of claim 1, further comprising permitting 
said at least one decay parameter to vary with at least one of 
said at least two generalized coordinates. 

4. The method of claim 1, further comprising: 

allowing an environment for said defining equations for 

said grid system to vary from a first time to a second 
time; and 

applying said method and said defining equations at the 
first time and at the second time to determine said 
solution of said grid system at each of the first and 
second times. 

5. A system for generating an elliptic grid of coordinates 
in two dimensions or in three dimensions, the system 
comprising a computer that is programmed: 


8 

to provide defining equations, valid near at least one 
boundary segment in a generalized coordinate system, 
of a selected grid system, where each of the defining 
equations has at least two independent Cartesian coor- 
5 dinate variables, has at least one generalized coordinate 
as a dependent variable, and comprises a partial dif- 
ferential equation, expressed in at least one generalized 
coordinate; 

] f) to provide a selected group of boundary constraints for the 

grid system, valid near the at least one boundary 
segment, where a decay parameter for at least one of 
the generalized coordinate dependent variables near the 
at least one boundary segment is determined as part of 
15 a solution for the grid system, rather than being pre- 
scribed initially; 

to provide defining equations and selected boundary con- 
ditions, having at least two independent coordinate 
variables and at least one dependent variable, for steady 
20 state heat transfer on a long thin fin, and to provide a 
correspondence between the at least two independent 
coordinate variables for the grid system near the at least 
one grid boundary segment with the at least two 
25 independent coordinate variables for the heat transfer 
problem; 

to provide a correspondence between a selected power of 
at least one heat transfer coefficient for the heat transfer 
problem and at least one decay parameter for the grid 
30 system near the at least one grid boundary segment; and 

to determine a solution of the grid system near the at least 
one grid boundary segment that incorporates at least 
one boundary constraint comprising the at least one 
decay parameter determined for the grid system; and to 
1 display the solution of the grid system. 

6. The system of claim 5, wherein said grid system decay 
parameter is directly proportional to said heat transfer coef- 
ficient. 

40 7. The system of claim 5, wherein said at least one decay 

parameter varies with at least one of said at least two 
generalized coordinates. 

8. The system of claim 5, wherein said computer is further 
programmed: 

to allow an environment for said defining equations for 
said grid system to vary from a first time to a second 
time; and 

to apply said method and said defining equations at the 
50 first time and at the second time to determine said 
solution of said grid system at each of the first and 
second times. 



